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Stability of recently predicted cage-like carbon cluster C$ (prismane) 
against its transformations to structures lying at lower energies is studied 
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theoretically. The smallest energy barrier inhibiting prismane transfor- 
mation is shown to be 0.44 eV, in qualitative agreement with the value of 
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activation energy for prismane decay calculated by molecular dynamics 
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simulations. A rather high value of the barrier suggests the possibility 
of experimental observation of this cluster at ambient pressure. 
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1. Introduction 

Since fullerenes Cgo have been discovered JL] and synthesized in macroscopic quan- 
tities @ , there is a growing interest in carbon clusters and nanostructures @, H; Hi • 
This interest has both fundamental and applied aspects. On the one hand, small 
carbon clusters composed of N ~ 10 — 100 atoms differ from macroscopic samples 

x; 

(and even from nanoparticles composed of N ~ 10 4 — 10 5 atoms) in that there are 
no "inner" atoms in such clusters, i.e. all atoms constituting the cluster lie on its 
"surface". A principal impossibility to subdivide interatomic bonding into bulk and 
surface parts results in essentially new physics and makes small clusters extremely 
interesting objects for basic research || [7J. On the other hand, some unusual prop- 
erties of small carbon clusters can be used for practical purposes. For example, 
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metastable cubane CgH 8 is considered as a high-energy- density material that can 
store a considerable amount of energy M; there is a hope that doping of C 36 fullerene 
can result in its phase transition to a superconducting state with an anomalously 
high critical temperature T c enhanced over the alkali-intercalated Ceo compounds 
§, etc. 

In order the small clusters could form bulk structures by means of interclus- 
ter bonding, they should be three-dimensional. Three-dimensional carbon clusters 
(cages) are believed to be stable for N > 20 only [TO], |TT|, while for clusters with 



iV < 20 the stable structures are either one-dimensional (chains) or two-dimensional 
(rings) |7], |9|, 11 1 . However, for iV < 20, carbon clusters can exist in three-dimensional 



metastable states whose binding energies are higher than the binding energy of the 
stable cluster with the same iV (for example, the occurrence of a metastable C14 



cage has been reported in Ref. [10), based on the results of density functional cal- 
culations). 



Several cage-like metastable isomers Cs have been proposed in Refs. [T|, [13], [14 
using the Harris-functional approximation combined with a generalized-gradient ap- 
proximation ||12j| , the local spin density and gradient-corrected approximations to 



the exchange-correlation energy [ 13 ] as well as the trasferable tight-binding potential 



14) . As far as we know, eight-atom isomers have not been observed experimentally 
yet, they are the smallest three-dimensional carbon clusters found so far theoreti- 
cally. The key issue for a possibility of experimental observation of such clusters 
is their relative stability. Note, however, that the stability of isomers Cg have not 
been quantified in refs. [113, [TJJ, e.g., the transformations to structures lying at lower 



energies have not been studied. Meanwhile, the results of ah initio calculations are 
known to depend strongly on the way that electron correlations are treated [0, [D| , 
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and hence a conclusion about the cluster met ast ability may appear to be just an 
artifact of approximations used. 

Recently we have predicted the existence of an eight-atom metastable cage- 



like cluster Cs that has the shape of a six-atom triangular prism with two excess 
atoms above and below its bases. We gave this cluster the name "prismane". It 
is shown in Fig.l. The binding energy of prismane is 0.45 eV/atom higher than 
the binding energy of the stable one- dimensional eight-atom cluster shown in Fig.2. 
Nevertheless, molecular dynamics simulations gave evidence for a relatively high 
stability of prismane, the activation energy for prismane decay was estimated to 



be about 1 eV [fDH , suggesting that prismane lifetime is rather long and that this 
cluster may be observed experimentally. 

However, although the finite-temperature molecular dynamics simulations provide 
the most direct way for evaluation of the cluster lifetime, such calculations are ex- 
tremely time-consuming since the process of cluster decay is probabilistic in nature, 
and hence a huge statistics should be accumulated in order to draw a definite con- 
clusion about the values of the cluster lifetime and activation energy. Alternatively, 
useful complementary information concerning characteristics of the metastable state 
can be provided by calculations of the heights of energy barriers inhibiting sponta- 
neous transformation of the cluster to the lower-energy atomic configuration. Such 
calculations have been carried out, e.g., to study the atmospheric-pressure stability 



of energetic phases of carbon [15] and polymeric nitrogen |16] . It was the purpose of 
this work to find the minimum energy barrier separating the metastable prismane 
structure C§ shown in Fig.l from the stable chain structure in order to see if such a 
cluster can be observed experimentally at non-exotic conditions. 
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2. Computational details 

To calculate the binding energy of a cluster having an arbitrary atomic configu- 
ration, we have used a transferable tight-binding potential [jll, [nfl that had been 



proven to reproduce accurately the energy-versus-volume diagram of carbon poly- 
types and to give a good description of both small clusters and bulk structures of 



carbon JD], [IT]. On the one hand, this potential describes the structure and ener- 
getics of small carbon clusters quite well, the difference in bond lengths and binding 
energies between our results and available ab initio calculations |7j usually did not 
exceed 10%. On the other hand, the technique used greatly simplifies a reliable eval- 
uation of characteristics of metastable states as compared with ab initio approaches 



H II- 



The search for the minimum energy barrier separating two locally stable atomic 
configurations is equivalent to the search for a saddle point in the total-energy sur- 
face as a function of the atomic coordinates. Such a saddle point corresponds to the 
unstable configuration of atoms in the cluster. The algorithm of the search for a sad- 
dle point is based on a continuous transformation of the cluster along the directions 
of atomic motion in the softest oscillation mode of the cluster, i.e. the mode with 
the minimum oscillation frequency. This algorithm is analogous to the "soft-mode 



walking" procedure p(| in which the information regarding both gradients and Hes- 
sian (energy curvatures) is used to systematically proceed up the bottom of a valley 
(local minimum) to a transition (saddle) state. In uphill walks from the minimum- 
energy geometry of the cluster to the transition state, the energy is maximized in 
the soft-mode direction and minimized in all orthogonal directions, thereby forcing 



the walk to be in a stream bed [20|. Apart from the efficient search for a saddle 



point, such an algorithm gives the energy spectrum of mechanical oscillations of the 



cluster. The knowledge of this spectrum allows for a reliable discrimination between 
the metastable and long-living unstable states of the cluster (sometimes this cannot 
be done by means of molecular dynamics simulations). 

All calculations were carried out at pressure P = 0, i.e. in fact at ambient 



(atmospheric) pressure 1T6 



3. Results and discussion 

We define the binding energy E b of a C 8 cluster as 

E b = E(8)-8E(l), (1) 

i.e. as the difference between the total energy of the cluster, E{8), and the energy 
8E(1) of eight isolated carbon atoms. The binding energy of the metastable prismane 
(P) is E^ = —40.80 eV. The binding energy of the stable eight-atom chain (C) is 
E£ = -44.37 eV. 

Starting with the prismane structure, we first computed the energy spectrum of 
its mechanical oscillations. We have found that all oscillation frequencies are real, 
giving evidence that the prismane structure corresponds to a local minimum in the 
total-energy surface as a function of the atomic coordinates, i.e. the prismane is 
indeed a metastable cluster. 

We now consider the energetics associated with the prismane transformation to 
the stable chain configuration in the soft-mode direction. Fig. 3 shows the difference 
AE b = E b — E^ between the binding energy of the transformed C 8 cluster, E b , and 
E^ versus the generalized transformation coordinate R defined as a mean square 
deviation of atomic coordinates X{, y,, Z{ of the transformed cluster from their cor- 
responding values Xio, yio, z i0 in the prismane structure (i — 1 — 8 is the number of 
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atom in the cluster): 

8 

R 2 = J2K x i ~ x io) 2 + (Vi ~ Viof + (zi - z i0 ) 2 ]. (2) 

i=l 

One can see that the function AEb(R) has a local minimum at R = corresponding 
to the metastable prismane structure. The value of AEb(R) increases with R and 
passes through the maximum SI at R 0.5 A. This maximum corresponds to 
the unstable configuration "scorpion" shown in Fig. 4. The binding energy of the 
"scorpion" is E^ 1 = —40.36 eV, so that the energy barrier inhibiting spontaneous 
prismane transformation to the lower-energy atomic configuration appears to be 
U = El 1 - E[ = 0.44 eV. 

We have found (see Fig. 3) that further increase in R leads to the prismane trans- 
formation at R ~ 6.5 A not into the stable eight-atom chain (as one could expect) 
but into another metastable configuration "frying pan" (FP) shown in Fig. 5. This 
configuration corresponds to the local minimum in the total-energy surface. The 
binding energy of the "frying pan" E FP = —42.85 eV is 2.05 eV below the binding 
energy E p of the prismane but 1.52 eV higher than that of the stable chain, E^ . 

Next we have studied the transformation path of the "frying pan" in its soft- 
mode direction. Fig. 6 shows the value of AEb = Ej, — E^ as a function of the 
transformation coordinate R defined by Eq.(||), where x i0 , y i0 , z i0 are now the atomic 
coordinates of the "frying pan". One can see that AEb(R) first increases with R, 
passes through the maximum S2 at R ~ 2.7 A and then decreases down to zero 
at R ~ 6 A, so that the stable chain structure (see Fig.2) is finally reached. The 
maximum of AEb(R) corresponds to the unstable configuration shown in Fig. 7. The 
binding energy of this configuration is E^ 2 = —42.08 eV, so that the energy barrier 
inhibiting the "frying pan" transformation to the linear chain is equal to U\ = 
E^ 2 — E^ p = 0.77 eV. Hence, the prismane transformation into the stable eight- 
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atom chain is a two-step process that takes place via an "intermediary" metastable 
configuration " frying pan" . 

Note that there also exist the paths of prismane transformation other than that 
shown in Figs. 3 and 6. The latter, however, is characterized by the smallest energy 
barrier U inhibiting transformation of the prismane structure to the lower-energy 
configurations. It is the path with the smallest barrier that is key to the issue of 
metastability. 

We stress that the data of finite-temperature molecular-dynamics simulation JT4 



are consistent with the results presented above. For all temperatures studied, 800 K 
< T < 1800 K, the prismane first transforms to the "frying pan" configuration, and 
next to the linear chain. The activation energy E a = 0.82 eV for prismane decay 
estimated in ]TJ| on the basis of the temperature dependence of prismane lifetime 
is about twice the height of the smallest energy barrier U = 0.44 eV calculated in 
this work. Note, however, that the cluster lifetime depends, in general, on the full 
topology of the total-energy surface, beyond just the minimum energy barrier. 

4. Conclusions 

We have numerically examined the metastability of recently predicted cage-like 
carbon cluster prismane Cg. It was demonstrated that the minimum energy bar- 
rier separating the prismane structure from the lower-energy atomic configuration 
amounts to 0.44 eV. Combined with the data on finite-temperature molecular- 
dynamics simulation of prismane decay, our result points to rather long lifetime 
of the prismane in its metastable state, suggesting the possibility of experimental 
synthesis of this cluster. 
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Figure captions 

Fig.l. Metastable prismane C 8 . Binding energy Eg = —40.80 eV. Bond lengths: 
d AB = 2.31 A, d AC = 1.28 A, d AD = 1.47 A. Bond angles: LB AC = 90°, LADB = 
104°, LABE = 60°. 

Fig.2. Stable chain C 8 . Binding energy Eg = -44.37 eV. 

Fig. 3. Calculated difference A£& = Eb — Eg between the binding energy Ef, of the 
cluster transformed from the prismane structure (P) in the soft-mode direction 
and the binding energy Eg of the stable eight-atom chain as a function of the 
transformation coordinate R. The prismane structure corresponds to R — 0. The 
saddle point configuration "scorpion" (SI) is realized at R m 0.5 A. The lower-energy 
metastable configuration "frying pan" (FP) is reached at R « 6.5 A. The energy 
barrier for the prismane — > "frying pan" transformation is U — Eg 1 — Eg = 0.44 
eV. 

Fig. 4. Unstable configuration "scorpion" (SI in Fig. 3). Binding energy E^ 1 = 
-40.36 eV. 

Fig. 5. Metastable configuration "frying pan" (FP in Figs. 3 and 6). Binding 
energy E[ p = -42.85 eV. 

Fig. 6. Calculated difference AEt> = Eb — Eg between the binding energy of 
the C 8 cluster transformed from the "frying pan" structure (FP) in the soft-mode 
direction and the binding energy Eg of the stable eight-atom chain as a function of 
the transformation coordinate R. The "frying pan" structure corresponds to R — 0. 
The saddle point configuration (S2) is realized at R fh 2.7 A. The stable chain 
configuration (C) is reached at R ~ 6 A. The energy barrier for the "frying pan" — > 
chain transformation is U 1 = Eg 2 - Eg p = 0.77 eV. 

Fig. 7. Unstable configuration S2 (see Fig.6). Binding energy Eg 2 = —42.08 eV. 
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